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Abstract 

We study the optical conductivity of a doped graphene when a sublattice symmetry breaking 
is occurred in the presence of the electron-phonon interaction. Our study is based on the Kubo 
formula that is established upon the retarded self-energy. We report new features of both the real 
and imaginary parts of the quasiparticle self-energy in the presence of a gap opening. We find an 
analytical expression for the renormalized Fermi velocity of massive Dirac Fermions over broad 
ranges of electron densities, gap values and the electron-phonon coupling constants. Finally we 
conclude that the inclusion of the renormalized Fermi energy and the band gap effects are indeed 
crucial to get reasonable feature for the optical conductivity. 
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1. INTRODUCTION 



There is a considerable interest in understanding the effects on properties of particle 
due to the interactions with environment, for instance the coupling of electrons to lattice 
vibrations or electron-phonon coupling. The electron-phonon coupling plays an essential 
role in the theory of high temperature superconductivity and they exist in other material 
such as nanotubes, Ceo molecules and other fullerenes pj. Also it is important to consider 
the electron-phonon coupling in transport properties. 

Graphene, a single layer of carbon atoms, [2J is disputable the first true two-dimensional 
lattices. Graphene is thermodynamically stable and there is indeed ripple structures on 
graphene sheets. Lattice displacements due to the ripple structures are symmetric with 
respect to their close carbon atoms and couple to the carrier densities. The electrons moving 
through the sheet are coupled to the out-of-plane phonons and therefore the electron-phonon 
coupling Pl ays an important role in the transport properties g fl Q. The coupling of 
electrons to out-of-plane optical phonons can be modeled by a Holstein type coupling [6J. 
In this model the coupling of electrons to dispersionless optical phonons is essentially local. 
The electron-phonon coupling has been carefully examined and has been shown to give rise 



to Kohn anomalies in the phonon dispersion at edge points in the Brillouin zone where the 
phonons can be studied by Raman spectroscopy 7, Q, Q]. An alternative strategy for the 
electron-phonon coupling measurement is based on the analysis of the G-peak linewidths 
and its broadening. 

The optical conductivity is one of the most useful tools to investigate the basic prop- 
erties of materials. Both the excitation spectrum of materials such gaps, phonons and 
interband transitions and the scattering mechanisms leave their distinct traces in trans- 
port. It was shown that the infrared conductivity of graphene is basically independent 



of the frequency [10. 
effect of e 



|l3| and experimentally confirmed this manner 14|, |15|. The 
ectron-phonon interaction in gapless graphene has been discussed by several au- 
thors u&HQQSi 

directed towards understanding this effect on the optical conduc- 
tivity. 

The energy spectrum of the Dirac electrons in a graphene layer that epitaxially grown on 
a SiC substrate has been measured by Zhou et al. 21] and they observed an energy gap of 
about 200 meV opened up in the electronic spectrum. They attributed the opening up of the 
gap is due to the breaking of the A and B sublattices symmetry 22| . The optical response of 
a gapped graphene is of important for an understanding of optoelectronic devices. Moreover, 
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the optical spectroscopy can be used for measurements of the magnitude of the energy gap. 

In this paper we consider the sublattice symmetry breaking mechanism for a gap opening 
in a pristine doped graphene sheet and study the impact of the electron-phonon coupling 
on the electronic conductivity of the electron-doped gapped graphene using Kubo formula 
at zero-temperature. We show that the renormalized velocity is suppressed due to the 
electron-phonon interaction. There is a shift in the chemical potential and we show that the 
interacting chemical potential is less than the noninteracting one due to the electron-phonon 
coupling. The optical conductivity is affected by Pauli blocking below twice value of the 
renormalized interacting chemical potential and gap values. 



2. MODEL HAMILTON! AN AND THEORY 

We consider the simplest form of Hamiltonian that describes the interaction of electron 
with an optical phonon mode, called the Holstein model. The honeycomb lattice can be 
consider in terms of two triangular sublattices A and B. We consider electrons in 7r-orbital 
of carbon atoms by using the tight-binding Hamiltonian in addition to the effect of the 
electron-phonon coupling due to localized Holstein phonons and a gap opening procedure 



due to sublattice symmetry breaking 24] . The total Hamiltonian in momentum space can 
be expressed as 



H = -tY^[4>(k)al a b kt(T + h.c] 

k,u 

+ D Yl XoWp, a a P +Kv + 6p,A+k^K c k + c -k) 

p,k,a 

k k 
k 

where ak CT or b^u is the fermion annihilation operator in k— space on sublattice A or B, 
respectively and t is the nearest neighbor hopping parameter [25j. The band gap, 2 A 
has a nonzero value as a result of breaks the symmetry between sublattices, A and B. 
We consider that the noninteracting chemical potential, fio be larger than the gap value 
representing the electron-doped system. The electron-phonon coupling is determined by 
D and furthermore u n denotes the fermionic Matsubara frequency. Moreover, Ck is the 
annihilation phonon operator, ujq is the frequency of the out of plane vibrations of the 
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optical phonon and Xo = y 2MNlu with M is ion's mass and N denotes the number of unit 
cells. 4>{k) = Y2s exp _l<5k with 5 being the vectors connecting the three nearest neighbors on 
the honeycomb lattice 26|. <p(k) reduces to hvpk/t in the Dirac cone approximation 25]. 

The matrix element of noninteracting Green's function with the gap of the electronic 
spectrum is determined by following expression 

AT„fl \ 1 



£ A (A;) / iho n + no- A£ A(fc) ' 

in which a,/3 = A, B and we have defined parameters Taa = —Tbb = A, Tab = T* BA = 
—t<f)(k). The quasiparticle excitation energy is f A (fc) = ^t 2 \(fi(k)\ 2 + A 2 . Note that at 
zero-temperature fiQ = £ a (&f) with kp is the Fermi momentum of charge carriers. 

An exact evaluation of the self-energy is only possible in some special cases. The matrix 
elements of the self-energy calculated to the lowest order in the electron-phonon interaction 
and is defined as 

ZaA^p) = -k B TY^D 2 xlD^\k,iv)G^Jp-q,iuj n -iv) (2) 



k.i' 



where G [ l a nd are the zero-order electron and phonon Green's functions, respec- 



tively 2y, |27|. In Holstein phonons, xo an d D^°\k,iu) are momentum independent and 



thus the phonon propagator is simplified by 

D<®(k,iu) = -2lu /(v 2 + u 2 ). (3) 

We restrict our calculations to the lowest order self-energy that is sufficient if Migdal's 
theorem, states that vertex corrections in the electron-phonon interaction can be neglected 
if the typical phonon frequencies are sufficiently smaller than the electronic energy scale, is 
valid. Therefore, we can neglect the vertex corrections since the self-energy is fc-independent. 
Using the contour integration, we can perform the summation over the bosonic frequency in 
the expression of the self-energy and finally the self-energy yields as 



fc,A=±l 



C A (fe) 



iVo + n F (\£ A{k) - no) | Np + 1- n F (\£ A{k) - /i ) ? ^ 



ihw n + hw - A£ A (fc) + fJL ihu n - ftw ~ A£ A (fc) + A*Q 

where JV = l/(e fao/,:BT - 1) and np(x) denotes the Fermi-Dirac distribution function. To 
calculate Ess(^n), the gap value A might be replaced by —A in Eq. HI It should be noted 
that Y^ABiiUn) = ^BAii^n) = in the Dirac cone approximation. The explicit expression of 
the self-energy will be computed in the following. 
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2.1. Finite doping with a gap opening 



We consider the low excited electron energy where the noninteracting electron spectrum 



energy is given by a/ [hv-pk) 2 + A 2 24| . To evaluate the zero-temperature retarded self- 
energy evaluated at the Fermi surface, we integrate Eq. H] over k and then decompose the 
results into E j (c<j) = Sq(^) + AE-'(w) where 

3te££(w) = (0W O ) { — £ In — — — - — — — + — In } 5 

Z7r vp (nw + /i ) — [rk^o + A) z vp oj- — u 

3mEj( W ) = - 7l ^-(gco Q ) 2 {-^^e(-u J+ - cu )e(tku + fi + hj + Jh 2 v 2 F k 2 + A 2 ) 

.Z7T 



+ Uj 2 ^ 6(cu_ - cu )e(-^ - /i + + v/ h 2 vpk 2 + A 2 )} 
f F 

here a^- = a; + (/xo + j A) / H with j = +1(— 1) refers to sublattice A(B). k c is the ultraviolet 
cut-off momentum |16j and finally the coupling constant g = \/NDxo/^o being the order of 
unity. The area of the unit cell is A c = a 2 3\/3/2 with a = 1.42A. The extra terms take the 
following form as 

UeAE j (u) = ^(^ ) 2 {-^ln| ( " + ^ ( "-"" 0) J - ^ In | / ~ ^ 2 |} (6) 
V 1 2vr vy u; L vp \u-u Q ){u_ + u y vp \ujj) 2 -uj 2]i w 

SmAS-'fu;) = — 7r^— ^(gu;o) 2 { — % — - 6(u;_ + co )Q(— to — cu ) 
2n vp 

_ ^r_^ 0( _ w + ^ (c_ - Wo )} (7) 

V F 

If A = 0, the self-energy reduces to massless Dirac graphene which addressed in Ref |l6| . 
Therefore, we have generalized the retarded self-energy expression to gapped graphene. 
Once the retarded self-energy is obtained, the quasiparticle properties of system due to the 
interaction of the electron-phonon can be calculated. The renormalized electronic spectrum 
is given by the Dyson equation as = ^(k)+dteTi(E] s ). Notice that according to the Dyson 
equation, we might distinguish the noninteracting chemical potential from the chemical 
potential of the interacting system due to the fact that 9fte£(u;) is not vanished for doped 
graphene when uj tends to zero. We thus have 

H = hq + 3fJeE(w)|^ = o . (8) 



The renormalized velocity, on the other hand, is given by 
v* Hv F k/U(k) + (v F H)- 1 d^eE(k, uj) 



\u)=0, k=kp 



(9) 



within the Dyson scheme 26] . The self-energy is independent of the momentum, accordingly 



its fc-derivative is zero. Consequently, the renormalized velocity is obtained analytically 

1 | ( 9Uo Y[]n\ ( A2 + (^ C ) 2 ) i 



v*(l + (A/hv F k F y/ 2 ) \v F k c J 1 ((ji + fko ) 2 - A 2 ) 

- (aio + jA + ^o) ^ + ^ q)2 _ A2 + 2^—}- (10) 

2.2. Optical Conductivity 

The optical conductivity can be calculated from the Kubo formalism. To this end, we need 
to obtain the current operator which is a composition of the paramagnetic and diamagnetic 
terms, i.e. j a = + j^p^p- We do need to modify the hopping parameter in the presence 
of an electromagnetic field [3] and then expand it up to the second order in the vector 
potential A(t). The current operator expressions do not change in the presence of the gap 
value and therefore by assuming that the electric field is in the direction of x-axis, we have 

f x = -i£ - 3)4 (k)b a (k) - (cf>*(k) - 3K(*)&t (*)] (11) 

<j,k 

where £ = tea/h and then the Kubo formula for conductivity is given by 

a fa) = 5 % > I A ^(^ + l jH (12) 
xx{ ' iA s (uj + ir]) ihA s (uj + ir]) 1 ' 



where A s is the area of sample and A xx (iu n ) = j^ kBT dre luJ " T < T T j^(r)j^(0) > 26]. 
We have ignored vertex corrections in the Kubo formula since we worked in nearly highly 
electron doped graphene for which the Dirac cone approximation is applicable. It was shown 
that the vertex corrections is essential for the low density carriers of the DC conductivity of 



graphene. 23j After a lengthy but straightforward algebra, we find 

A f kc f°° d e 
5mA«(w) =( 2 -^ / kdk —(n F (e + u)-n F (e)) 

x {{2t 2 \ ( f ) {k)\ /l )A AB {k,e)A AB {k,e + u;) 

+ (9-|0(A;)| 2 ) [A AA (k,e)A BB (k,e + u) 

+ A BB (k,e)A AA (k,e + u;)}} (13) 

where the spectral functions are the imaginary part of Green's function which take the 
following forms: 



Here 



Q± = ifiw n + /i ± A, ^AA(BB) = ^+(-) ~ ^BB(AA)(^n) 

and <&ab — *&ba — 1- The integral over k in Eq. [13] can be performed analytically and 
accordingly one dimensional integral will be needed to be calculated numerically Note 
that the interacting chemical potential is used instead of the noninteracting one because 
of the nonzero value of £j(0) . It should be noted that by setting A = 0, the optical 
conductivity results are different with the results given in Ref. [l6] due to the fact that we 
have implemented the interacting Fermi energy in the formalism. 



3. NUMERICAL RESULTS 

We have considered the system with the phonon energy being huo = 0.2 eV [jjj]. Although 
the order of coupling constant is unity, we consider a larger value to seek its effect better. 
We have found that the value of the quasiparticle properties for sublattices A and B are 
different at most about 0.8% due to the gap opening. We will then present only the results 
of the sublattice A. 

In Fig. 1, we have shown the results of the real and imaginary parts of the retarded 
self-energy for the electron-doped system, (/x > A) at n = 5 x 10 12 cm~ 2 . SmE(w) vanishes 
in \u\ < ujq at which point it jumps up to a finite value because only then can a quasiparticle 
decay by boson emission. It drops towards the zero for uj < —ujq and then increases linearly 
showing a marginal type physics which happens in the Coulomb electron interactions in 
undoped graphene 29]. Notice that SmS is not symmetric with respect to change of the 
sign of frequency. In addition, 9mS(a;) vanishes when \ fvuj + fru} + fi \ < A due to the effect 
of the gap opening and 9mS tends to zero at — huo — /Jo for A = 0. These behaviors can 
be determined explicitly from expressions given by Eqs. Eland [7J 

In Fig. lb we can see logarithmic type singularities 



at uj = ±u) Q and uj = —uj 



(/io — A)/H for the results of 3fte£. The extra singular behavior is due to the gap effect. 
It should be noted that the singularity at u — ±cj would be washed out if a momentum 
dependence of phonon spectra is used. In addition, there is a cancelation of the logarithmic 
singularity at u — —ujq — (/io + A)/h. The logarithmic singularity can be determined to 
the argument of the logarithm in Eqs. E] and We have obtained an expression for the 
interacting density of states too through the spectral function. The singularities manner 
lead to kink structures in the interacting electronic density of states. In the results, there 
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are three kink structures in the interacting density of states where one of them is associated 
to the gap. The kink structures would affect to physical quantities and transport properties 
through the interacting electronic density of states. 

The renormalized velocity as function of the densities, gap values and the coupling con- 
stants are shown in Fig. 2. The renormalized velocity is suppressed due to the electron- 
phonon interaction and the gap values too. We have found a nonmonotonic behavior of 
v* with respect to the electron density when the gap value increases and results are shown 
in Fig. 2b. At small gap values, v* decreases with increasing density however it changes 
behavior at large gap values and behaves like conventional two-dimensional electron sys- 
tems. Therefore, we expect that the electron-phonon interaction renormalized the electronic 
quantities at the Fermi surface by a factor v*/vp. 

The optical conductivity scaled by o"o = e 2 /4h as a function of energy for different values 
of (a) the coupling constants and (b) the gap values are shown in Fig. 3. First of all, a 
tends to a minimum value at ujq. Moreover it basically increases around uj > ujq due to the 
contribution of the Holstein phonon sideband. In the case of noninteracting electron-phonon 
system, a has a sharp structure, step function manner, at 2/xq due to interband transitions 



and the conductivity increases by a factor o 
frequencies decreases and approaches to a 



two, a = 2<7q at u = 2 A and finally at higher 



13J. By switching interaction on, the chemical 



potential becomes weaker and consequently the position of the sharp structure changes to 
2/j, which is smaller than 2/xo- This behavior is clearly shown in the Fig. 3 which did not 
consider in results discussed in Ref. 16| . At g = 0, the conductivity is larger than ctq about 



2/x and then tends to <To in gapped graphene. However, a always remains smaller than <jq 
in gapless graphene. The gap dependence on the optical conductivity is shown in Fig. 3b. 
First, the gap opening makes the chemical potential bigger therefore the sharp structure in 
the a tends to larger u values. Second, the scattering mechanism increases by increasing 
the electron densities and then the optical conductivity changes and becomes smaller. 

Another point of interest for experiments is the density dependence ( in units of 10 12 
cm -2 ) of the optical conductivity ( Fig. 4) as a function of frequency at 2A = 0.2 eV. Note 
that the noninteracting chemical potential values associated to the electron densities used 
in Fig. 4 are /i = 0.154,0.279,0.382 and 0.831 eV, respectively with giving A = 0.1 eV. 
The optical conductivity increases by increasing the electron density around u however a 
decreases faster by increasing the density at high frequency. The sharp structure of the 
optical conductance tends to higher frequency by increasing the electron density. The sharp 



8 



position occurs at 2/x which is always smaller than 2/xq for the same system. 



4. CONCLUSION 



we have calculated the optical conductivity of gapped graphene, including the effect of the 
lowest order self-energy diagram due to the electron-phonon interaction by Holstein Hamilto- 
nian. We have reported an extra logarithmic singular behavior associated to gap value in the 
real part of the self-energy We have found the density, gap value and the electron-phonon 
coupling dependence of the renormalized velocity and the interacting chemical potential. 
The optical conductivity is affected by these physical quantities and Pauli blocking below 
twice value of the renormalized chemical potential and the gap values. We conclude that 
the inclusion of the renormalized Fermi energy and the band gap affects are indeed crucial 
to get reasonable feature for the optical conductivity. The gap dependence of the optical 
conductivity would be verified by experiments. 
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FIG. 1: (Color online) Imaginary (a) and real (b) parts of the self-energy as a function of energy 
evaluated at Fermi energy for different gap values at the coupling constant g = 3.0 and density 
u = 5.0 x 10 12 cm~ 2 . 




FIG. 2: (Color online) Renormalized electron velocity for (a) the different values of coupling 
constants at n = 5.0 x 10 12 cm~ 2 , (b) the different value of density (in units of 10 12 cm -2 ) at 
g = 3.0. 
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FIG. 3: (Color online) Optical conductivity as a function of energy for (a) the different values 
of coupling constant at A = 0.1 eV and (b) the different value of A at g = 3.0. We consider 
n = 1 x 10 13 cm~ 2 . 
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FIG. 4: (Color online) Optical conductivity as a function of energy for different values of density 
( in unites of 10 12 cm" 2 ) at g = 3.0 and A = 0.1 eV 
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